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Abstract. In a frame of quasi-crystal approximation the dispersion equations are obtained for 
the wave vector of a coherent electromagnetic wave propagating in a media which contains a 
random set of parallel cylinders with possible overlapping. The results are compared with that 
for the case when a regularity at the cylinder placement exists. 

Introduction 

Under certain conditions the electrochemical anodization of metals leads to the formation of 
the oxide layer with parallel cylindrical pores on the surface of metal. The radius of the 
cylinders is much less than the thickness of the layer. These oxide films are used as 
corrosion-protective, strengthening, decorative coatings and attract much attention in 
microelectronics, in particular, for the creation of non-linear and active thin-film elements, 
and also for substrates of large hybrid microcircuits. Porous oxides are perspective for 
multilevel systems of interconnections. In dependence on the conditions of anodization and 
kind of metal it is possible to obtain films with different degree of ordering in pore placement 
of pores. The placement of pores is maximally chaotic when the pores are placed randomly 
and could be overlapped. The overlapping pores form the cavities of complex configurations. 
The placement of pores is maximally ordered, when the points of the cylinder axes 
intersection with the plane perpendicular to the axes lay in the nodes of some lattice, for 
example a hexagonal one. An investigation of the optical properties of oxide films allows to 
find a degree of pores ordering. If the wavelength is much larger than the radius of pores, the 
optical properties of the porous layer can be described by the effective refractive index. This 
index is to be found from the dispersion equations. 

The dispersion equation for the effective refractive index of the system of parallel non- 
overlapping cylindrical fibres in media was deduced in the work [1]. The equation allows to 
describe the optical properties of the layer with the regular pores placement. In present work 
we obtain the dispersion equation for a medium with chaotic distribution of pores (the pores 
are located randomly and can be overlapped). We also analyse an influence of regularity on 
the effective refractive index. 

1. Dispersion equation for a medium with completely chaotic cylinder 

placement 

Let us consider a system of parallel dielectric cylinders. We assume that if the distance 
between the centres of cylinders is less than distance 2R (R being the radius of the cylinders), 
they form the common scatterer with the constant refractive index. If the concentration is 
small, the number of overlapping cylinders is also small. So the majority of scatterers are 
cylindrical. If the concentration of cylinders increases, the number of the overlapped 
cylinders increases. However, even for medium consisting of such complicated scatterers it is 
possible to derive the dispersion equation for a wave vector using the statistical averaging 
technique [2-4]. Without the loss of generality we can set to unity the dielectric constant of 



the surrounding media and the dielectric constant of the scatterers set to £ . Let the z-axis is 
directed along the cylinder axis, then, the positions of cylinder axis at the plane perpendicular 
to z are specified by two-dimensional random vectors s ...s N . The dielectric constant £ of 
N+l overlapping cylinders can be written in the form 



£(p,z)=£(p) = ! + (£-!> 
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where 0(p) = < . The expression (1) is done so, that £ (/>) = £ , if there are one or 
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more cylinders in the given point. The electromagnetic field satisfies Maxwell's equations, 
whose solutions can be written in the form: 

E(r,s ,s v ..s N ) = E(p,s ,s v ..s N )e lk ' zZ . 

For E(p,s ,s 1 ...s N ) we obtain the equation: 

DxDxE(p,s ,s 1 ...s N )-k 2 E(p,s ,s 1 ...s N ) = (£-l)k 2 £(p)E(p,s l ...s N ), (2) 
where D = (d /dp,ik' z ). Let average this equation over positions of all the scatterers except 
s , in point of view that the cylinders are independent and distributed with the density of 
probability f(s t ). So we multiply (2) by f(s l )...f(s N ) and integrate by ds v ...ds N . At 
\p — S 1 < R we obtain: 

DxDx(E(p,s Q ))-k 2 (E(p,s Q ))=(£-l)k 2 (E(p,s )). (3) 
At \p — S \> R the averaging of right part (2) gives: 
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where £l(p) = $ f(s)G(p-s)d 2 s . 



Here by analogy with the quasi-crystal approximation [5,6] we have put 

(E(p,s ,s l ...s N _ l ,s N )) - (E(p,s ,s l ...s N _ l )) - ... - (£(/»,«!>). 



Let us go to a limit of infinite medium by directing N — > °° under Nf(s) — n , 
where n is the number of scatterers on unit area, then 

N£l = n 7tR 2 . (5) 

Taking into account the uniformity of averaged medium we write down [5]: 

E(p,s)) = e ik ^(E(p-s,0)). (6) 



A vector k' = (k' L ,k' y ) is the desired wave vector [4], connected with 
effective refractive index by n eff =k'/k . Substituting (5) and (6) into (4) and uniting with (3) 
we obtain: 

DxDx{E(p,0))-k 2 (E(p,0))=(£-l)k 2 e(p)(E(p,0)) + 

+ (£-l)k 2 (l-0(p))n' o e ik ' ±p §e~ ik ' lS (E(s,0))d 2 s , (7) 
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n' = (l-e ' n ^ tRl )l%R 1 . (8) 

The value n' tzR is the filling parameter, i.e. a fraction of an area in a cross section 
occupied by cylinders. Equation (7) allows to find the dispersion equation for a wave vector 
k' . At complete filling of a medium by cylinders n — > °° , n' — > 1/jtR . It is easy to see 
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that in this case the solution of (7) is (E( p,0)) = E exp( ik\p ), k\ +k' z = ek . Let 
rewrite (7) in integral form, using the Green function [7] 
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where D = \ — , ,ik' \,V = Jk 2 -A:' 2 , ReV>0, H L and J L are the Hankel 
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function of the first kind and the Bessel function, respectively. 
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As seen from the equation (7) at \p \< R , (E(p,0)j satisfies to the wave equation 

D x D x (E (p,0)) - ek 2 (E (p,0)) = . 
Solution of this equation can be represented as [7,8] 

(E(p,0)) = Ua T :E™ap) + a T m E E:*ap)}e im<p , (11) 



where A = -yjek — k' z . The expressions for E™ and E TE are given in [1]. 

Substituting (11) into (10) just as in [1], and projecting (10) on the vectors 

e^+k^m/ Rv 2 e z and e z , multiplying by exp(-/m<p), and integrating by d(p at 
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\p\ = R , we obtain an infinite-dimensional system of equations relative of a n and a n : 
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where [HJvRjJJ k\R)]= k\RHJvR )J'J k\R )-vRH'JvR)JJ k\R ) and so 

on. Everywhere in (12) we imply summation over subscript n. 

The condition of solvability of the system (12) gives the dispersion equations for wave 
vector k' , which are found simpler than equations for non-overlapping cylinders obtained in 
[1], since it is possible to bypass a calculation difficulty with of the infinite-dimensional 
determinant of the system. The system (12) has two types of the solutions symmetry. The first 
one corresponds to an ordinary wave. 
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Taking into account the symmetry (13), we obtain the dispersion equation in the form 
which is suitable for the numerical calculations by iteration method and Newton method. 
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For an extraordinary wave 
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we have the dispersion equation 
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where D n = a 2 n - ji J n . 

In static limit kR — > from the equations (14) and (15) it follows: 
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where £ n = l + n' 7tR 2 ( £ -1 ), £ ± =l + 2(e -l)n' 7lR 2 /((£ + l)-n' 7lR 2 (£-l)) are the 

dielectric constants of the medium with cylinders embedded for the transverse and 
longitudinal electrostatic field correspondingly [8]. 

Let us show how to apply the dispersion equations in the case of porous media. 
Parameters of the problem are the values of the wave vector in the medium k m and inside a 

cylinder k c . Instead of k m and k c we used k = k m and £ = (k c / k m ) . For porous 

medium, k m = J£ m k and k c = k , where £ m is the dielectric constant of the material, and 



k = CO I c is the wave vector in vacuum. From this, it follows that k = - s J£ m k and 
£ = 1 / £ m . The effective refractive index of the porous medium is given by the expression 
n e ff=k'/k , where k' is found from the dispersion equation solved with parameters 



k = J£ m k and £ = l/£ m . 
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Let us consider the case of oblique incidence of an electromagnetic wave on the layer 
of parallel cylinders, which axes are perpendicular to the layer surface. For this case 
k' ± = k ± , where k ± is the component of a wave vector of incident wave, which is 

perpendicular to the z axis (parallel to the layer surface) and the relation k ± I k is sine of the 

angle of incidence. At given k' ± , the k' z has to be obtained from the equation (14), if the 

wave is polarized perpendicular to the plane of incidence, or from (15) if the field strength 



vector lays in the plane of incidence. Let us show the results of numerical calculations for the 
porous medium with n m =U6. The dependence of the real and imaginary parts of the 

effective refractive index n eff = \k\+k' z f / k on the normalized frequency (i.e. on the 
parameter COR I c ) is shown in Fig. 1,2 for the incident angles 0°, 60° and for different 
porosity parameter 7tn' R . At COR I c — > the real part of the effective refractive index goes 
to the static value. The imaginary part has a maximum at k R ~ 1 . 

2. Comparing with the case of the regular cylinder placement 

System of non-overlapping cylindrical pores already possesses some degree of ordering. This 
degree of ordering increases if the concentration of pores becomes larger. It is possible to 
increase degree of ordering without changes in the concentration in the following way: let us 
imagine that each pore is surrounded by impenetrable cylindrical shell with the radius R s . For 

completely ordered media with hexagonal cells R s =l/^3v2n' . For the function w(x) in 
the equation (20) of work [1] we should take the following function 

R. 
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where w a ' (x) is a radial distribution function for the rigid disks on a plane with the filling 

parameter a' = 7M' R s . 

In Fig. 3 the frequency dependence of the effective refractive index is plotted for 
different degrees of ordering of pores. The wave is propagating in direction parallel to the 
axes of cylindrical pores. The curve number 1 corresponds to pores inside the non-overlapping 



shells with the radius R s =0,81/-y/3V2n' . The porosity parameter 7Tn' R for all curves is 

the same. Only the degree of ordering is changed. If the ordering becomes more chaotic, the 
imaginary part of the effective refractive index caused by the attenuation of the coherent wave 
becomes larger. It is maximal for the completely random pore placement with possible 
overlapping. The real part of the refractive index is also strongly dependent from ordering. 

Conclusion 
We obtained in analytical form the dispersion equations (14 and (15) for the wave vector of 
the electromagnetic wave propagating in a media which consists of randomly placed parallel 
cylinders with possible overlapping. In this media the maximally chaotic placement of the 
scatterers is realised. Together with the obtained in [1] equation (20) for the non-overlapping 
cylinders our equations allow to find the effective refractive index for any degree of medium 
regularity. The imaginary part of the refractive index increases when the placement of pores 
becomes more chaotic. 
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Figures 
Fig. la, b. Dependence of the real (a) and imaginary (b) parts of the porous medium 

o 

(n m =l,76) effective refractive index on the normalized frequency for the incident angle 

( parallel to the cylinder axes) , the porosity parameter equals: 1 - 0.1; 2 - 0.2; 3 - 0.3; 4 - 0.4; 
5-0.5. 

Fig. 2a, b. Dependence of the real (a) and imaginary (b) parts of the porous medium 
( n m =1,76) effective refractive index on the normalized frequency for the incident angle 60° , 

the porosity parameter equals: 1 - 0.1; 2 - 0.2; 3 - 0.3; 4 - 0.4; 5 - 0.5. The solid curve 
corresponds on the extraordinary wave, the dashed one corresponds on the ordinary wave. 

Fig. 3a, b. Dependence of the real (a) and imaginary (b) parts of the porous medium 
(n m =\&5) effective refractive index on the normalized frequency effective, the porosity 

parameter equals to 0,2 . The curve (1) corresponds to the pores, placed in non-overlapping 
circles with radius R S =1,72R , (2)- non-overlapping «hard» pores, (3)- chaotic placement of 

o 

pores with possible overlapping. Light incidence angle . 
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